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ABSTRACT 



An interactive computer model of a highly enriched pressurized water 
reactor was developed, using the applicable plant parameters from the 
Shippingport Atomic Power Station. The point reactor kinetics equations 
for one delayed neutron precursor group were linearized using small 
perturbation theory. The model included both moderator and Xenon-135 
reactivity feedback effects, as well as an automatic reactor protection 
and average reactor coolant temperature control system. The thermal 
response of the model plant was simulated for normal operating transients 
induced either by control rod or turbine load changes. The post shutdown 
Xenon transient response was also modeled. The interactive program was 
coded in FORTRAN-IV language, and the simulation program was coded in 
IBM CSMP-III language. 
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I. INTRODUCTION 



Computers have been used in the design and analysis of nuclear 
reactors since the inception of reactor technology, and many codes have 
been developed. More recently, digital computer programs have been used 
as learning devices for nuclear engineering students. The formulation 
of the computational problem of predicting the behavior of a nuclear re- 
actor has been greatly facilitated with the use of high level computer 
1 anguages . 

The purpose of this work was to develop a computer assisted learning 
device to be used by students taking nuclear engineering courses at the 
Naval Postgraduate School. The program graphically displays the simulated 
kinetic and thermal transient responses of a pressurized water reactor 
power plant. The reactivity feedback effects of Xenon-135 poisoning and 
moderator temperature are incorporated into the model. Normal operating 
transients, starting from a steady state critical condition, can be 
induced by ramp changes in either control rod position or turbine load. 

The initial reactor power level and reactivity change mechanism are 
chosen by the program user. For user convenience, these inputs are 
prompted and entered interactively, after which the simulation is run 
with no further user action required. 

A pressurized water reactor (PWR) was modeled as it is the most 
common type of reactor used for power plant application. As a realistic 
reference, the model design incorporated the applicable characteristics 
of the Shippingport Atomic Power Station, whose nuclear parameters were 
the most compatible, of available PWR core data, with the assumption of 
a highly enriched core made in the model design. 
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The model was developed by treating the reactor kinetics and other 
plant component thermodynamic relationships as transfer functions [1]. 
The point reactor kinetics equations were linearized using a Taylor ex- 
pansion or small perturbation technique [2]. This same perturbation 
technique and lumped parameter analysis were applied to the plant's 
linear differential heat transfer and Xenon-135 equations [3]. The 
individual transfer functions were then integrated into an overall plant 
block diagram. The program also features a reactor protection and 
average reactor coolant temperature control scheme. 

IBM CSMP-III language was chosen to formulate the model simulation 
as it has inherent routines to invert the transfer functions to the time 
doma i n , 
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II. MODEL DESIGN CONSIDERATIONS 



Lumped parameter analysis and first order perturbation theory were 
used throughout the model development. All the nuclear variables which 
were used in the reactor kinetics and Xenon-135 decay equations were 
considered to be averaged values of the variables over the neutron 
energy spectrum. Similarly, the thermodynamic variables which were 
used in the plant's heat transfer equations were considered to be 
averaged values over the volume of the particular component. 

Because of these approximations and other assumptions made in the 
model development, the simulation should not be considered a design or 
stability analysis. Instead, the simulation shows the model's large 
scale reactor kinetic and thermal trends during normal operating 
transients . 

For similar reasons, although certain plant parameters of the 
Shippingport Atomic Power Station were incorporated into the model, the 
simulation cannot be considered to reflect the operating characteristics 
of this power plant. The Shippingport core is a seed and blanket type 
whereas the model core is a uniform mixture of highly enriched fuel 
and support materials. 

A. PRESSURIZED WATER REACTOR 

The pressurized water reactor (PWR) is the most widely used reactor 
type in central power plant applications and the only type currently 
used in naval propulsion. 

A simplified schematic of the modeled PWR plant is shown in Figure 
1. The modeled reactor contains a highly enriched core which is light 
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Figure 1. Schematic of the Model Plant 

water moderated and cooled. The modeled plant contains two closed loop 
thermodynamic systems coupled by a heat exchanger. 

In the primary system, heat generated from thermal fission is trans- 
ferred to the coolant as it passes through the reactor, raising the 
coolant's temperature. The high temperature coolant leaving the reactor 
is circulated by a pump through tubes inside of a heat exchanger where 
it gives up some heat to a secondary fluid. A pressurizer maintains 
the primary coolant at a sufficiently high pressure so that the bulk 
temperature of the coolant is kept below the saturation temperature. 

On the secondary side of the heat exchanger, the temperature of the 
entering feedwater is raised to the boiling point and saturated steam 
is produced. The steam is delivered to a turbine, is condensed and 
returned to the heat exchanger, completing the second closed loop. 
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B. MODEL AND PROGRAM FLEXIBILITY 



The model was designed to simulate only normal operating transients 
from an initial steady state condition. Reactor accidents and startup 
were not considered in the model development. However, the model's re- 
actor control module will simulate either a full or constant insertion 
of control rods if certain parameters are exceeded. This feature was 
not incorporated for accident analysis, but to keep parameters within 
the limits of the assumptions made in the model development. 

The applicable Shippingport plant characteristics are fixed in the 
simulation program. User inputs are limited to choosing initial power 
level, from which other initial parameters are adjusted, and the plant 
perturbation mechanism, either control rod movement or turbine load 
change. These perturbations occur at fixed rates which limit the amount 
of reactivity that can be inserted during the simulation. 

Originally it was envisioned that an overall program would feature 
not only interactive input capability but also produce real time graphi- 
cal displays of the transient for a user at a time sharing computer 
terminal which had graphics capability. However, the computer language 
(IBM CSMP-III) required to solve the model's algorithms was not avail- 
able on the time sharing system, and furthermore, the programs long 
execution time precluded any real time response. 

In order to still provide as much user facility as possible with 
these restrictions, the interactive feature was partially retained for 
data input. Using a computer language available on the time sharing 
system (FORTRAN-IV) , a program which interactively prompts the user to 
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enter specified controlling parameters was developed. This program in- 
corporates logic routines which permit inputs only within requested 
ranges. This feature allows for data reentry if a user error is made, 
and ensures that only inputs compatible with the model development's 
assumptions are used in the simulation. 

After the input is completed, the initial transient conditions are 
displayed at the terminal. A separate internal control program then 
transfers the user supplied inputs into the simulation program. The 
control program then transfers the simulation program to the batch 
processing system for execution. The interactive program notifies the 
user that this has been done. Hard copy plots of the transient are 
subsequently produced. The post-shutdown Xenon behavior module is 
located in the interactive program and produces real time graphical 
displays at the terminal. 

Thus there are three distinct programs: 

1. An interactive program to prompt and receive user input, 
and also simulate post-shutdown Xenon behavior. 

2. A batch processed program which simulates the transient 
response to the user inputs. 

3. An internal control program which interfaces these two 
programs. 
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III. MODEL DEVELOPMENT 



A. POINT REACTOR KINETICS EQUATIONS 

A lengthly and formal derivation of the point kinetics equations is 
found in Reference 4 and will not be repeated here. More generally these 
equations are obtained from one group neutron diffusion theory with the 
assumption that the neutron flux is deparable in time and space, and 
with the inclusion of delayed neutrons [5], They are listed below. 

0(t) + ^A.,c,(t) t q(t) (1) 



dC.(t) 

~dt 




A .C.(t) 



i = 1 , 2 ,... 



( 2 ) 



where 


0(t) 


= 


Flux amplitude function 




/}{t) 


= 


Reactivity 




c.(t) 




Effective concentration of the ith delayed 
neutron precursor group 




Q(t) 


= 


Extraneous delayed neutron source strength 






— 


Effective delayed neutron fraction from 
the ith group 




^i 


= 


Decay constant of the ith group 


and 




= 


21 = Total effective delayed neutron 
i fraction 



These equations are referred to as the point reactor kinetics equa- 
tions not because the reactor is considered as a single point in their 
application, but since spatial variations are neglected, and 

are assumed to be constant. 

If the neutron flux spatial variation is assumed to be time invari- 
ant, 0(t) can be considered to represent the number of neutrons in the 
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core. Equation (1) may then be seen as a neutron rate equation where 

the three terms on the right hand side represent the rate of production 

of prompt, delayed, and source neutrons respectively. 

Equation (2) is a rate equation for the ith delayed neutron precursor 

group. The two terms on the right hand side represent production and 

decay rates respectively. The precursors are comprised of approximately 

thirty isotopes which historically have been divided into six groups 

-1 235 

with decay constants ranging from 0.0124 to 3.0 seconds for U. In 
the model, the precursors were considered to be represented by one 
effective group characterized by an averaged decay constant [2]. 




and an effective delayed neutron fraction 

Z/5 

In equations (1) and (2), reactor power may be substituted for neutron 
flux if the precursor concentration is modified by C = E^ ^ f^old 

since 

where P = Reactor power 

E^ = Energy released per fission 

= Macroscopic fission cross section 

0 = Integrated one group neutron flux 

For a reactor operating at power, the source term contribution to 
the overall neutron population is negligible, i.e., Q » 0. 
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Incorporating these assumptions, equations (1) and (2) reduce to; 



dPM = + Ac(t) 



( 3 ) 



dC(t) ^ 
dt 



/3 



P(t) - \C(t) 



(4) 



1 . Zero Power Reactor Transfer Function 

Consider an initially critical steady state reactor at some power 
level P^ at t 5 0. Equations (3) and (4) each yield: 

C P 

0 AA 0 



(5) 



Now let P(t) ■ Pq <^(t) 

c(t) = Cq + y*c{t) 
/3{t) + <^(t) 



( 6 ) 



where the zero subscript denotes the initial steady state value 
and the delta prefix a smal 1 perturbation imposed at t = 0 about 
this value. 

Noting that = 0 for a critical reactor and substituting equa- 
tions (5) and (6) into equations (3) and (4) yields for t ^ 0 



^ d'Ht) = <CoM - p(t) + X/c(t) 

^ <rc(t) = (/’pct) - 

where the term in equation (7) has been neglected. 

Taking the Laplace transforms of equations (7) and (8): 

s ^(s) = 0^(5) - -^</*P(S) +. Ac/C(S> 

n 

</’C(s) s - A</*C<5) 

A 



(7) 

( 8 ) 



(9) 

( 10 ) 
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Solving equation (10) for C(s), substituting this into equation 
(9), and rearranging, yields: 



c/*PCs) „ 



S+A 



■) 



( 11 ) 



This is the zero power reactor transfer function, so called because 
the reactor is assumed to be operating at a sufficiently low enough 
power that no feedback effects are realized. These effects are examined 
in the following sections. 

Equation (11) may be represented in block diagram form as shown in 
Figure 2. 



tCols) 






cTP(S) 

> ■ 



Figure 2. Block Diagram of the Zero Power Reactor 
Transfer Function 



where Z(s) = 



s (a 






S+X 



-) 



2. Limitations of the Point Reactor Kinetics Equations 

The point reactor kinetics equations derivation is based on 
the assumption that the spatial dependence of the neutron flux is 
negligible. This assumption limits the validity of these equations to 
transients where this remains a reasonable approximation such as those 
which result in only small changes in reactivity. The small perturba- 
tion technique used in the development of the zero power reactor 
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transfer function also requires only small changes in reactivity if the 
first order approximation is to hold. Specifically, must be less 
than 0.5/3 . Physically, when />is greater than /3 the reactor is 
critical on prompt neutons alone. The simulated transients imposed on 
the model were limited to ensure that excessive amounts of reactivity 
were not introduced. 

B, REACTIVITY FEEDBACK MECHANISMS 

In the zero power reactor point reactor, the power level is assumed 
to be so low that it does not affect the reactivity, thus there are no 
feedback effects. However, for a reactor operating at a useful power 
level, feedback effects do exist and the reactivity becomes an implicit 
function of the reactor power level (or neutron flux). This dependence 
arises since reactivity depends on macroscopic cross sections which in- 
volve the atomic number densities of material in the core: 

2 = No- 
where 2L = Macroscopic cross section (cm“^) 

N = Atomic number density (atoms/cm^ ) 

O' = Microscopic cross section (cm ^ ) 

The atomic number density can depend upon the reactor power level 
since the concentrations of certain nuclei are constantly changing due 
to neutron interactions. Material densities also depend upon tempera- 
ture which is a function of reactor power level and hence the flux.^ 



^Duderstadt, J. J. and Hamilton, L. J., Nuclear Reactor Analysis , 
1st ed., Wiley, 1976. 
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The two feedback mechanisms considered in the model were Xenon-let and 
moderator temperature. 

Reactivity can also be changed directly by an external source such 
as control rods containing a neutron absorbing material. Thus the over- 
all reactivity can be written as; 



where 



P 

^T= Fj P 

= Overall core reactivity 
= Externally added reactivity 
<^x = Xenon-135 feedback reactivity 

= Moderator temperature feedback reactivity 
Fx = Xeonon transfer function 
Fj = Moderator temperature transfer function 



(12) 



Equation (12) with the previously developed zero power reactor 



transfer function, given by equation (11), are incorporated in the 



block diagram below; 




Figure 3. Block Diagram of the Reactor Transfer Function 
with Xenon-135 and Moderator Temperature 
Feedback Loops, and External Reactivity 
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The next step in the model development was to derive the Xenon-135 
and moderator temperature transfer functions. 

1 . Xenon-135 Transfer Function 

Xenon-135 is the most significant fission product poidon because 
of its enormous thermal neutron absorption cross section and relatively 
large fission yield. This isotope is not only produced directly from 
fission but also from the decay of other fission products as shown in 
Figure 4. 




Since the decay of Iodine-135 and Xenon-135 (^^^Xe), with 

the largest half-lives, are the controlling steps in this decay scheme, 
it was simplified by making the following assumptions: 

1) All I is produced directly from fission (the production of 

Antimony-135 (^^^Sb) and subsequent decay to is considered 

instantaneous) . 

2) The short lived metastable Xe is ignored. 

1 35 

3) The removal of I by neutron absorption is negligible for 

the neuton flux levels used in the model (10^^ neutrons/cm^sec) . 
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4) All the decays to ^^^Xe. 

With these assumptions the effective decay scheme is shown in Figure 



5 . 




Cs 



Xe 



cr*0 

a 



Figure 5. Simplified Xe-135 Decay Scheme 



135 



1 35 

Using this decay scheme, the resulting rate equations for I and 
Xe are: 

(13) 

VxZfe5(t) + X(t) (14) 

135 3 

where I = I number density (atoms/cm ) 



1 35 3 

X = Xe number density (atoms/cm ) 



Zf = Macroscopic fission cross section (cm" ) 

2 

0 = Integrated one group neutron flux (neutrons/cm sec) 
Xj = Effective fission yield 

= Effective ^Xe fission yield 

135 1 

decay constant (sec ) 

135 “1 

= Xe decay constant (sec" ) 

X ,35 

^ = Xe microscopic thermal neutron absorption cross 

^ section (cm 2 ) 
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Again, using the first order perturbation technique, let 
I(t) = + cPlit) 

X(t) = Xq + </X(t) (15) 

0(t) = 13^ + cPUt) 

where the zero subscript denotes an initial steady state value and 
the delta prefix denotes a smal 1 perturbation about this value. 

Upon substituting equations (15) into equations (13) and (14), and 

taking the Laplace transform, the relationship between the perturbation 
1 35 

in Xe and 0 is derived (See Appendix A): 



<A(S) +M(yr?f X«) 



(f0<.s^ S^+ (oa*^otXx<-Aj)s+ Ax) 

The change in reactivity caused by a small perturbation in "^''Xe 
concentration is also derived in Appendix A: 



= Gx Cs) ( 16 ) 
135, 



where 



(/X(s) 



Xo + ^ 
135, 



= = 



U 



(17) 



X = Equilibrium Xe number density before the 
° perturbation (atoms/cm^) 

o oc 3 

U = Uranium-235 ( U) number density (atoms/cm ) 



= Microscopic thermal neutron absorption cross 
section of U (cm ) 

A perturbation in neuton flux ( c/*0) is directly proportional to a 
perturbation in reactor power ( c/*P) as shown by: 
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where = Energy released per fission 

o n r T 

= Macroscopic fission cross section of (cm’ ) 
or 

(/gics) _ 1 

cy'p(s) “ 

The product of equations (17), (16), and (18) yields: 



( 18 ) 



g^x(s) (fXLS) 0^0(5) _ g^(S) 

(/‘XCS) <f0(s) c/'P(S) </*P(s) 






(19) 



Equation (19) is the transfer function relating reactor power and 
135 

Xe feedback reactivity. This equation is shown in block diagram 
form in Figure 6 where it has been incorporated with the previously de- 
rived zero power transfer function given by equation (11). 




Figure 6. Block Diagram of the Reactor Transfer 
Function with Xenon-135 Feedback and 
External Reactivity 



2. Moderator Temperature Feedback 

Reactivity feedback from changes in moderator temperature occurs 
as a result of changes with moderator density. The density is also a 



24 



function of pressure, however, the pressure coefficient of reactivity is 
typically two orders of magnitude smaller than the temperature coeffi- 
cient, and was therefore not considered in the model development. The 
moderator density affects the moderator number density and hence the 
macroscopic scattering cross section of the moderator. 

The primary mechanism for the thermal ization of the prompt and de- 
layed neutrons is by elastic scattering interactions with the moderator 
nuclei. Hence, variations in the macroscopic scattering cross section 
will affect the rate at which neutrons become thermal i zed. Changes in 
the thermal ization rate affect the fission rate, or reactor power level. 
Power level changes affect the moderator temperature, thus a feedback 
loop is created. 

Because the moderator is also the coolant in the model, its temper- 
ature is not only a function of reactor power but also a function of the 
heat transfer process occurring in the heat exchanger, thus complicating 
the feedback loop. 

Therefore, in order to develop this feedback mechanism analytical ly, 
it is first necessary to model the plant's heat transfer processes. 

This thermal analysis is done in the following section. 

Another temperature feedback mechanism is the broadening of the 
Uranium-238 resonance absorption cross section for neutrons with increas- 
ing temperature. Because a highly enriched core was assumed in the 
model, with little Uranium-238, this effect was not considered. 
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C. THERMAL ANALYSIS 



1 . Reactor Heat Transfer Function 

A lumped parameter model was assumed. This simplification pro- 
vided a set of ordinary differential equation which were sufficiently 
accurate for the simulated normal operating transients. In the lumped 
parameter model, heat transfer in the reactor was assumed to occur at a 
single point. Thus, spatial variations were neglected. The core was 
considered to be a homogenized mixture of the uranium alloy fuel, the 
fuel cladding, and other structural materials, with a constant thermal 
capacity. 

The equation for the heat flow from the core was obtained from a 
basic heat balance. The heat generated from fission equals the heat 
required to change the temperature of the core materials plus the heat 
transferred to the coolant. On a per unit time basis, the heat balance 
is : 



where P(t) = Total Power generated in the core (Btu/sec) 

Tp(t)= Average temperature of the core materials (°F) 

TAv(t) = Average reactor coolant temperature (°F) 

Cp = Total thermal capacity of the core materials (Btu/°F) 

hpn^ = Total heat transfer coefficient (Btu/° F sec) 

A similar heat transfer balance must also hold for the heat being 
transferred to the coolant and transported out of the core. This heat 
is the last term of equation (20) and is transferred to the coolant as: 
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l^F«['^FCt:)-TAv(t)] = Ch + mmC [Th,(0 (21) 



where 



C^. = Total thermal capacity of coolant in the core 
^ (Btu/°F) 




= Coolant mass flow rate (Ibm/sec) 

= Average reactor coolant outlet temperature (°F) 

= Average reactor coolant inlet temperature (°F) 

= Specific heat of reactor coolant (Btu/lbrn^F) 



For simplification, the average reactor coolant temperature is 
assumed to be given by; 



Tav (t) = [ Tc • (jk) + Tho (t)] / 2 
Rearranging equation (20) yields: 



where 



TpC-tl + 't'l 



7" _ Cf 

•l “ TT“ 
•^FM 



d Tp^-tl 
dt 



TavU) + ii 

Cp 



P(t) 



( 22 ) 



(23) 



Substituting equation (22) into equation (21) to eliminate yields: 



TavW = 



where 


to - 








mnC 




11 


Cm 




^FM 



dt 

and 



TpCtl + 2 -^7^; (fc) 

•o 1 



( 24 ) 
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Now let 



P(t) = + 


d'P(t) 


Tp(t)= Tp + 


d'TM 


0 


Av'‘)= * 

0 




H Too * 




0 0 


'■'o 


‘Ap (t) 



( 25 ) 



where the zero subscript denotes a steady state value and the 
prefix a small perturbation about this value. 

Substituting equations (25) into equations (22), (23), and (24), 
eliminating the average temperatures Tp(t) and T^y(t), and taking Laplace 
transforms, the following relationship between reactor coolant inlet and 
outlet temperature perturbations is derived (see Appendix B); 



c^Tho(s1 = (26) 

where 

^ 



This equation is the reactor heat transfer function for the reactor 
coolant outlet temperature as a function of the reactor coolant inlet 
temperature and reactor power. 

A block diagram representation of this transfer function is shown 
in Figure 7. 
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/PW 




Figure 7. Block Diagram of the Reactor Heat 
Transfer Function 



In Figure 7 



Gc(S) = 






G^^(S) = 






Z. Heat Exchanger Transfer Function 

As in the derivation of the reactor heat transfer function, a 
lumped parameter model was assumed for the heat exchanger. Two points 
of energy storage were assumed, the primary coolant water and the water 
on the secondary side of the heat exchanger. The thermal capacity of 
the heat exchanger metal was included with that for the secondary water, 
since the thermal resistance on the primary side is predominant. A 
further assumption was made that the time spent by the primary coolant 
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while it passed through the heat exchanger was negligible in comparison 
with the time spent in the primary piping. 

The state of the steam produced on the secondary side of the heat 
exchanger was assumed to always be dry and saturated and that the 
secondary water is always at the saturation temperature for the existing 
pressure. These assumptions were justified because moisture separators 
and recirculation can provide high quality steam and preheating of the 
feedwater. 

With these assumptions, the following equations were obtained from 
a heat balance per unit time: 

= c„ ^ hT„[TAve(t)-Ts(4.)] (27) 



hT«[TAVg(tl -Ts(tl ] = Cs 



(28) 



where 



THi(t) 

Tco<‘> 

.V 

‘8 

Ts(‘) 

•'tm 

S 

P|_(t) 



m 



m 



Heat exchanger coolant inlet temperature (°F) 
Heat exchanger coolant outlet temperature (°F) 
Heat exchanger average coolant temperature (°F) 

Saturated steam temperature (°F) 

Total heat transfer coefficient (BTU/°F.sec) 

Total thermal capacity of coolant in the heat 
exchanger (BTU/°F) 

Total thermal capacity of heat exchanger metal 
and secondary water and steam (Btu/°F) 

Power delivered by the heat exchanger (BTU/sec) 
Coolant mass flow rate (Ibm/sec) 
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The total heat transfer coefficient, a function of the primary coolant 
flow rate (a constant in the model) and the heat transfer characteristics 
of the heat exchanger, was assumed to be constant. Also, the time delay 
in transferring heat across the heat exchanger tubes was neglected. 

This changed the shape of the initial thermal transient but had little 
effect on the basic dynamics of the secondary loop. 

Again, as a simplification, the average temperature of the coolant 
in the heat exchanger was assumed to be given by 




I 



+ ) 



(29) 



The power delivered by the steam generator is proportional to the 
product of the steam flow rate and the difference in enthalpy between 
the steam and feedwater 



= m3 (Hs“Wfw) +■ Cs 

where m^ = Steam flow rate (Ibm/sec) 

= Saturated steam enthalpy (Btu/lbm) 
Hp = Feedwater enthalpy (Btu/lbm) 



This equation assumes the steam and feedwater flow rates are always 
equal and thus neglects any instabilities in the secondary steam. 

As discussed in Reference 4, the enthalpy of saturated steam is 
nearly constant over a wide range of pressure, varying from 1198 Btu/lbm 
to 1204 Btu/lbm over a pressure range from 200 to 800 psig. The feed- 
water enthal py, which depends on condenser pressure, is usually between 
50 and 100 Btu/lbm. Thus the enthalpy difference may be regarded as a 
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constant, and the power delivered by the heat exchanger is directly 
proportional to the steam flow rate. 

The impedance to steam flow caused by the turbine is nearly indepen- 
dent of turbine speed. If constant backpressure is assumed, the steam 
flow rate is directly proportional to the throttle opening at a given 
pressure. Thus, 



where 






^s = Saturated steam pressure (Ibf/in'^) 

A = Proportional ity factor which is a function 
of the throttle setting 1 bm i n^ 

Ibf sec 



The numerous assumptions made in the heat exchanger model development 
limit the accuracy of the resulting equations. However, the errors in- 
volved in these assumptions are usually less than the amount of uncertainty 
in the engineering valves of the coefficients used in the equations. 

Recalling the previous assumption that 

—it — ' ° 

Substituting equation (29) into equation (27), yields: 



TcoCt)] = Tco(t)] - Ts(t)J (30) 

Solving equation (30) for T^ , gives: 






2Ts(fe) - Th^ (1-Ki) 



( 31 ) 
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where 



Z m^c 
Ktm 



Substituting equation (29) into equation (28) gives: 

Ktm j'i [rH:iW+ -Ts(0 . = Cs + Pt(t) 



Rearranging this expression. 



Cs cJTs(fc) ^ ^ |-[TH.(t) + Tc«(o] + -^ Pc(t} 



HtM dt 



(32) 



Let 

(t) = + cf'T^it) (33) 

0 °o 

Ts(t) = + c/'T3(t) 

PL(t) = Pl + c/'p^_(t) 

^0 

where the lowest zero subscript denotes a steady state value and 
the delta prefix a smal 1 perturbation about this value. 

By substituting equations (33) into equations (31) and (32), and 
taking the Laplace transform, the following set of equations is derived 
(see Appendix C) : 



</’rs (s) = 



F [ Cs) t (s) ] - (s) 

Ts + 1 



(34) 






2</'’fs (s) - (5) [l - Kll 

1+ Ki 



( 35 ) 



where 
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These equations are represented in a block diagram in Figure 







Figure 8. Block Diagram of the Heat Exchanger 
Transfer Function 

With c/’Tu and c/’P, as inputs, equations (34) and (35) form an 
Hi L 

algebraic loop in </Tp and </t<.. These two equations were solved 

Lo i 

with the use of an inherent functional routine available in the CSMP-III 
language which was used to formulate the model. 
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3. Primary Piping Transfer Functions 



While circulating through the primary loop, the coolant under- 
goes mixing and transport delay effects. For example, a transient in 
the coolant temperature at the heat exchanger outlet does not appear at 
the reactor inlet until sometime later. Where there are volume or flow 
direction changes, as in the reactor coolant inlet plenum, mixing occurs 
causing a smearing of a temperature transient. Both of these effects 
were approximated in the model by combinations of two types of time 
delays: a pure transport delay and a simple mixing delay. 

Assuming no mixing in the primary piping and no heat loss with 
perfectly insulated pipes, pure transport delays are encountered in the 
piping runs between the reactor and heat exchanger. With these assump- 
tions, the temperatures involved in the transfer of heat between the re- 
actor and heat exchanger can be given as 



(t) = (t - T 3 ) 

where the inlet coolant temperature to the heat exchanger T^. has the 

same form as the outlet temperature of the reactor after a fixed 

0 

transport delay This transport delay can be approximated by the 
following differential equation which is derived in Appendix D. 



Th. * "^3 

1 



dTHi 

dt 



= T, 



(36) 



Similarly, the transport delay from the heat exchanger outlet to 
the reactor inlet plenum is*. 



Tc. (t) + 

ip 



dT. (t) 

ip 

dt 



= T. (t) 
^0 



(37) 
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where 



(t) = Reactor inlet plenum coolant temperature (°F) 

ip 

Tp (t) = Heat exchanger outlet coolant temperature (°F) 

^0 

This can be approximated by 

cJTc.p(t) ^ 

Tc^pCt) 1^ — ^ — = TcoCt) (37) 



As in the previous development, let 



T 



Ho 



(t) 




(t) 



Tc (t) 



Tp (t) 

^ip 




(t) 

(/t (t) 

^•p 

/Tp (t) 

0 

Ac (t) 

• ^ip 



(38) 



where again the lowest zero subscript denotes a steady state value 
and the delta prefix a small perturbation about this value. 

After substituting equations (38) into equations (36) and (37), 
taking the Laplace transform, the following transfer functions are 
derived (see in Appendix D). 



(S) _ 1 



(39) 



(/Tc^p (s) _ 1 

c/’Tco(s) ~ H-HS 



(40) 
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Combined mixing and transport effects are encountered at the reactor 
and heat exchanger inlet and outlet coolant plenums. For simplification, 
only the mixing effects at the inlet plenums were considered in the 
model. Assuming perfect mixing and after performing a heat balance on 
the plenum concerned, the combined mixing and transport effects are ex- 
pressed by differential equations developed in Appendix D of the form 

+ T^(t) = T,(t) 

where 

T. (t) = Plenum inlet temperature (°F) 

Tg(t) = Plenum outlet temperature (°F) 

T = Mixing delay 

As before, by using a small perturbation technique, and taking the 
Laplace transform, the following transfer function is derived (see 
Appendix D): 



<Aho (S) 






(5) 1 +TS 

Block diagram representations of the coolant piping transport and 
mixing transfer functions is shown in Figure 9. 
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D. OVERALL PLANT BLOCK DIAGRAM 

The transfer functions developed previously for the zero power point 
1 35 

reactor, Xe feedback, reactor and heat exchanger heat transfer, and 
the primary piping were interconnected resulting in the overall model 
block diagram shown in Figure 10. The moderator temperature feedback 
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/P(s) 




^Thp(s> 





f 






[ 

/ 










Figure 9. Block Diagram of the Transport 
and Mixing Delay Functions 



loop is seen to consist of the heat transfer process and primary piping 

transfer functions which yield and </Tm . These temperatures 

^i ^0 

are summed, then halved, yielding T^^, which is then multiplied by 
the negative temperature coefficient ®^T, generating the moderator 
temperature feedback reactivity /^j. 



E. MODEL CONSTANTS 

The following plant parameters from the Shippingport Atomic Power 
Station were obtained from Reference 6 and used in the model development. 
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Figure 10. Overall Plant Block Diagram 



1 . General Parameters 



a. Reactor thermal power 

b. Reactor coolant system pressure 

c. Reactor coolant average temperature 

d. Steam pressure at full load 

2 . Reactor Coolant System 

a. Reactor coolant flow rate 

b. Reactor coolant outlet temperature 
c« Reactor coolant inlet temperature 
d. Coolant volume in core 

3. Reactor Core 

6.8 ft 



0 

4. Nuclear Data 

a. Thermal neutron flux 

b. Prompt neutron lifetime 

40 



a. Configuration 

b. Size 

c. Fuel load (seed) 

d. Composition (seed) 

1 ) Water 

2) Fuel alloy 

3) Zircalloy 

e. Control Rods 

1 ) Total rod worth 

2) Scram time 



231 MW 
2000 psia 
523° F 
600 psia 

6280 Ibm/sec 
538° F 
508° F 
103 ft^ 

Right cylinder 
dia X 6 ft. high 
75 kg 

43.5 v/o 

30.0 v/o 

34.1 v/o 

0.256 k 

35 sec time delay 
1 .0 sec rod drop 

X 10^^ n/cm^ sec 
5.6 X 10'^ sec 



c. Effective delayed neutron fraction 0.0077 

d. Temperature coefficient of reactivity -3.1 X 10’*^ </k/°F 

5. Reactor Protection Setpoints 

a. Scram 

1) High reactor power 138% 

2) High reactor coolant outlet temperature 550°F 

b. Cutback 

1) High reactor power 114% 

2) High Startup rate 1.74 Decades per min. 
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IV. RESULTS 



A. INTERACTIVE PROGRAM 

Figures 10-16 are examples of the interactive program output. As 
shown in Figure 10, the user is first given a description of the purpose 
of the overall program and general instructions for input entry. As 
shown in Figure 11, the user is then prompted to enter an initial steady 
state power level, after which the corresponding major plant parameters 
are displayed. The user is then prompted to choose the type of simula- 
tion, either plant transient at power, or post-shutdown Xenon-135 
behavior. The user must be at the Tektronix 4012 console in the computer 
center if the latter is chosen. 

If the post-shutdown Xenon-135 behavior is chosen to be examined, 
graphs such as those shown in Figure 13 are generated and displayed on 
the Tektronix' s screen. In addition to showing the time response of the 
Xenon, the time of its peak and associated maximum reactivity are dis- 
played. 

If the plant transient simulation is chosen, the user is then prompted 
to choose the mechanism for initiating the transient, either control rod 
movement or turbine load change. In Figure 14, the user has chosen the 
former. The program then prompts the user to enter the direction and 
time length of the control rod movement. In Figure 15, the turbine load 
change was chosen and the program request the final turbine load. In 
either case, a summary of the transient inputs is displayed with a 
notice that the interactive portion is complete. 
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Figure 11. Interactive Program, Introduction and Instructions 
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Figure 12. Interactive Program, Initial Power Level and Simulation Choice Inputs 
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Figure 13. Post-shutdown Xenon-135 Behavior 
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Figure 14. Interactive Program, Control Rod Movement Inputs 
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Figure 15. Interactive Program, Turbine Load Change Inputs 



Figure 16 is an example of the interactive program's logic routines 
which will only permit inputs within the requested ranges to be accepted 
for the simulation. 

B. SIMULATION PROGRAM 

The program generates Versatec plots showing the plant's power, re- 
activity, and temperature transient responses for the user's inputs. 
Figures 17-20 are composites of these plots. 

In Figure 17, the transient was initiated by a simulated 20 to 40 
percent ramp change in turbine load at 1/2 percent per second. In 
Figure 18, the transient was initiated by a simulated 60 to 40 percent 
change in turbine load in the same manner. Both of these figures show 
the characteristic power demand following and inherent stability response 
typical of a PWR plant with a constant average temperature program and 
negative temperature coefficient. 

The inherent stability feature is further displayed in Figures 19 
and 20. In Figure 19, the transient is initiated by a simulated 10 second 
inward movement of the control rods at a reactivity insertion rate of 
1.25 X 10’^ (Pk per second from an initial 50 percent power level. In 
Figure 20, the transient response to an outward control rod movement 
with the same parameters is shown. As seen in both cases, the reactor 
power stabilizes about the initial turbine load and the average reactor 
coolant temperature stabilizes at a new level to compensate, via the 
negative temperature coefficient, for the reactivity inserted by the 
control rods. 
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ENTER AN INITIAL STEADY STATE POWER LEVEL BETWEEN 10* AND 100* PERCENT. 
NOTE! A HINIMUN OF 10% POWER IS REQUIRED TO OPERATE THE REACTOR COOLANT 
PUMPS AND OTHER PLANT AUXILIARIES 
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Figure 16. Interactive Program, Logic Routine Example 
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Figure 17. Power, Reactivity, and Temperature Response 
to a 20 to 40% Turbine Load Change 
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Figure 18. Power, Reactivity, and Temperature Response 
to a 60 to 40% Turbine Load Change 
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Figure 19. Power, Reactivity, and Temperature Response 
to a 10 second Inward Control Rod Movement 
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Figure 20. Power, Reactivity, and Temperature Response 
to a 10 second Outward Control Rod Movement 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The model responds to normal operating transients with the charac- 
teristics of a PWR plant. The interactive program provides the user with 
the facility to initiate the simulation or examine a real time display 
of post-shutdown Xenon-135 buildup and decay. 

As previously discussed, the purpose of this work was to develop a 
learning device for nuclear engineering students at the Naval Postgrad- 
uate School. The result was a relatively simple model of a PWR power 
plant and a simulation program that suffers from a long execution time. 

The model's simplicity resulted mainly from the consideration of 
an effective single group of delayed neutron precursors in the reactor 
kinetic equations and the use of lumped parameter analysis in the plant's 
heat transfer processes. The model's sophistication could be increased 
by: 

1. Expanding the reactor kinetic equations to consider six 
delayed neutron precursor groups. 

2. Developing a multi-section heat transfer model for the 
reactor and the heat exchanger. 

3. Incorporating variable reactor coolant mass flow rate and 
heat transfer coefficients. 

These additional features will also increase the complexity of the 
simulation program, thus aggravating the long run time problem. This 
concern might be eased by taking a different approach to the model's 
formulation than the transfer function method. While this method was 
readily coded using the CSMP-III language, physical insight to the plant 
dynamic processes was lost when the Laplace transform and subsequent 
grouping of constant terms was performed. The use of state variable 
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theory might not only allow this physical appreciation to be retained 
but also result in shorter program run times. 

Regardless of the formulation method, a real time response to the 
simulation can probably only be achieved by using an analog computer. 
The wide range of time constants associated with the equations (prompt 

_5 

neutron life-times on the order of 10 seconds to Xenon decay half 
lives on the order of hours) require a small numerical integration time 
interval over a large time period. The result is a prohibitively long 
run time on a digital computer for an interactive program. 

It is hoped that refinement of the model will be continued. While 
the existing model does reflect the general characteristics of a PWR 
plant to normal operating transients, it does exhibit relatively large 
power over/under shoots in response to turbine load changes as seen in 
Figures 17 and 18. 
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APPENDIX A 



DEVELOPMENT OF THE XENON FEEDBACK TRANSFER FUNCTION COMPONENTS 



1 . Derivation of G (S) 

A 

135 135 

The previously developed rate equations for I and Xe from the 
1 35 

effective Xe decay scheme shown in Figure 5 are: 



dl 



^ =. - Ax ret) 



(A 1) 



ZfiS M + AiJ (t) - AxX(t) ^(t) XIO 



(A 2) 



135 3 

where I(t) = I number density (atoms/cm ) 

135 3 

X(t) = Xe number density (atoms/cm ) 

Sf = Macroscopic fission cross section of U (cm" ) 

0 = Average integrated one group flux (neutrons/cm sec) 

^I = Effective ^^^I fission yield 

^X = Effect! ve^^^Xe fission yield 

^I = ^^^I decay constant (sec~^ ) 

Xx = ^^^Xe decay constant (sec~^) 

135 2 

(j*X = Xe thermal neutron absorption cross section (cm ) 

I(t) 



Let 



Iq + </l(t) 



X(t) = + <A(t) 

0(t) = + (/0{t) 



(A 3) 



where the zero subscript denotes an initial steady state value and 
the delta prefix denotes a small perturbation about this value. 
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With the reactor at an initial steady state with a flux level 0^, the 

135 135 

equilibrium values of I and Xe are found by setting their time 

dependence in equations (A 1) and (A 2) equal to zero. Thus, 



Xo = 



Xo = 



Xf $25o 

XI 

Zf <^a t Ai I< 



Substituting equation (A 4) into equation (A 5) 



(A 4) 



(A 5) 



X. = 



Ax + Oi^^o 



(A 6) 



Substituting equations (A 3), (A 4), and (A 6) into equation (A 1) 
and (A 2), neglecting the cPx term (first order approximation) 



d. 

d-t 






Sf - Ajc/'X 



(A 7) 



i C/-X = t (^xZf-q 3 'Xoic/’jZ(-(A,+aaV«)</X (A 8 ) 

Taking the Laplace transform of equations (A 7) and (A 8) and solving 
for </l(S) and <fx(S) 



cPl(5) = 



S rhz 



</X(s) = 



Arc/^I (SI t C^x - Oa^Xq) cP^Cs) 

(s + Ax tCaVo) 



(A 9) 
(AlO) 
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Substituting equation (A 9) into equation (AlO) 



c/'XO) = 



-h (s»Ai)(yxS^ -c3a^Xo)o^C3) 
C5+Xx)( S+A.X+- 



Expanding, collecting terms, and solving for 



cy^xisi 

c/‘<2(s> 



c/’Xfs) S + Ax (^x +Vx2p - 05 ^X>) _ 

s — ^x^LS) 

s*+ (Oa^^o + Ax tAi)s ■*■ ^x) 

This is equation (16) on page 23. 



2 . Derivation of 

The effective core neutron multiplication factor k is given by 
the familiar "six factor formula" found in the literature 






where 



^ - Thermal fission factor of the fule 

6 = Fast fission factor of the fuel 

p = Resonance escape probability 
f = Thermal utilization factor in the core 
PpNL ■ ^3St neutron non-leakage probability 
PfNL “ Thermal neutron non-leakage probability 



(A 11) 



The core reactivity is defined as 

k-1 

^ ■ k (A 1.2) 

Consider a reactor at an initial critical steady state condition 
1 35 

with an equilibrium Xe concentration Xo. By the definition of 
critical ity 
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In this condition 



ko = 1 

y^O = 0 



(A 13) 



fo 





(A 14) 



where 




= Fuel macroscopic thermal neutron absorption cross 
section 



21 



core 

3o 



= Initial core macroscopic thermal neutron absorption 
cross section 



^Xo 

^3 



1 35 

= Initial Xe macroscopic thermal 
cross section 



neutron absorption 



and the macroscopic thermal neutron absorption cross sections of the 
moderator and other core materials has been neglected, a reasonable ap- 
proximation in a highly enriched core. 

1 35 

Now impose a small perturbation in the Xe concentration on the 

1 35 

core. Neglecting the effect of the Xe perturbation on neutron leakage, 
a reasonable assumption for a large reactor, 

X = Xo + </x 
k = ko + <^k 
f = fo + </f 

/> = /^o + = cf/^ since /><, = 0 

The remainder of the terms in equation (A 11) are unchanged, thus 



and 



cfV = <ff 



ko 




since ko = 1 
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where 



■P = 



z; 

zr^ 



si 

Sa“%05Vx 



from equation (A 12) and (A 13} 



(A 15) 



.-J 



(A 16) 



substituting equations (A 14) and (A 15) into equation (A 16) 



= 1 - 






Vx 



<^a 



Vx 



57 oot-< 



7COV-C 



oa'^xo+ sr 



or 



</*x 



OgX 



= 



This is equation (17) on page 23. Since fuel depletion effects are 
not considered during the short simulated transients, the fuel number 
density U is assumed constant and invariant in the time and the 

Laplace domains. 
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APPENDIX B 



DEVELOPMENT OF THE REACTOR HEAT TRANSFER FUNCTIONS 



The previously developed equations describing heat transfer from the 
reactor to the coolant are 



where 



and 



TfU) = TAv(t) + P(t) 



^ 2 . •%] t I 2 s Tp (t) + 2 4 ||Tc^ (t1 

lAvcti = 



(B 1) 

(B 2) 
(B 3) 



Tp(t) 



^Ho(^) 

Tci(t) 



Average core material temperature (°F) 

Average reactor coolant temperature (°F) 

Reactor coolant outlet temperature (°F) 

Reactor coolant inlet temperature (°F) 

P(t) = Total power generated in the core (Btu/sec) 
p = Total thermal capacity of core materials (Btu/^F) 

= Total thermal capacity of coolant in the core (Btu/°F) 
C = Specific heat of coolant (Btu/blm°F) 
m_ = Coolant mass flow rate (Ibm/sec) 



m 



^FM 


= Total 1 


io 


Cm 


m^C 


ri 


„ Cf 

hpn 


'(z 


= £h- 

Wfm 
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Let 



P(t) > p„ 


+ 


(/p(t) 




+ 


cy’lp(t) 


' ^AV 


+ 




-- ^Ho 

0 


+ 




Tc.(t) = Tc. 

1 10 


+ 





where the lowest zero subscript denotes a steady state value and 
the delta prefix a small perturbation about this value. 

From substitution of equations (B 4) into equations (B 1), (B 2), 
and (B 3), with the reactor in a steady state at t ^ 0, 

TFo * TaVo Po 

TAVo<i + 2 (B 5) 

TaVo = 4- (Thp^ -h TcjJ 

and 

(B6) 

(Af(o) - c^7av(o) = cAIho(o) => (o) = O 

Impose the perturbations at t=0. Substituting equations (B 4) and 
(B 5) into equations (B 1), (B 2), and (B 3) with the initial conditions 
given in equations (B 6), 

c/TpCt) + fl £ = ^TAvCt) •h^c/’PCt) (B 7) 

«/TAv(fcl [l + 2 ^ [</Tav C t)] = ^ ^ Ct) (B 8 ) 

c/*lAv(tl * % [^g/’Tho (t) + (tl^ (B 9) 
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Taking the Laplace transform of equations (B 7), 
with the initial conditions given by equations (B 6) 


(B 8), and (B 9) 


Tf( 5 ) [sfl+l] = 1av(5) -f- -^PCS) 


(B 10) 


Tav(s) [57*2 + 142 = Tf(s) + 2 Tc(s) 


(B 11) 


TaV(S) = [ Th (S) 4 Tc (5)j 


(B 12) 



where the delta prefix and lowest subscript have been deleted for 
readabil ity 

Substituting equation (B 12) into equations (B 10) and (B 11) 



Tf (sYltl) = ( Th + To) + i p (B 13) 

HTh + Tc) (stj+i + 4^) = Tp + aiTc (8 ’“I 

where the s domain dependence notation has been deleted for 
readability. 

Solving equation (B 13) for Tp and substituting into equation (B 14) 

j(Th + Tc) (sf2+n.2^)(sfin) = -|(THtTc) + i P + 2 ^(Sfl+l)Tc 



After expanding the products and collecting terms 

After multiplying through by and solving for T^( (s)) 



</Tho(S) 







|s-l 


}c/'Tc.(5)+ ^(Pp(s) 








[I4 




s+i"! 
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where 
This i 



-r 

= —. 3 3nd all notation has been restored. 
12 

equation (26) on page 28. 
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APPENDIX C 



DEVELOPMENT OF THE HEAT EXCHANGER TRANSFER FUNCTION 
The previously derived heat exchanger heat transfer equations are 






g~s(fe) - TH^(fc) [l- ><il 

[l + Ki] 



(c 1) 



^ + TsCt^ = T [ THi(t) + ICo^) ] + hTM (C 2) 

K-tm o't 

where exchanger coolant inlet (°F) 

T^^(t) = Heat exchanger coolant outlet (®F) 

T^(t) = Saturated steam pressure (PSIQ) 

Pl(T) = Power delivered by heat exchanger (Btu/sec) 

Cj = Total thermal capacity of heat exchanger 

metal and secondary water and steam (Btu/°F) 

C = Specific heat of coolant (Btu/lbm°F) 

hj,^^ = Total heat transfer coefficient (Btu/®F.sec) 

m^ = Coolant mass flow rate (Ibm/sec) 



and 



=: 



Zm-nC 

Kth 



Let 



0 

+ cfljCt) 



(C 3) 



Ts(t) =T3 
PL(t) = Pl 



C^PL(t) 



where the lowest zero subscript denotes a steady state value and the 
delta prefix a small perturbation about this value. 
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With the plant in an initial steady state at t < 0, equations (C 1) 
and (C 2) are 

2Tso - ±+ Kj. 



T<=Oo = 



and 



[l + ICi] 

Ts^ = I ( Th-^ fTco^) + ^ Plo 

cPTu-^(o) = cPTqq(,o) =. c/’Ts(o) = o 

Impose the perturbations substituting equations (C 3), (C 4), and 
(C 5) into equations (C 1) and (C 2) 



(C 4) 



(C 5) 



(C 6) 



<rrc = 



2</'Ts -ci%i 



(itKi) 



(c 7) 



£ ^ A(t) t c/%(0 = 1 Acw] - (C 

Taking the Laplace transforms of equations (C 7) and (C 8) with the 
initial conditions given by equation (C 6) 

I [<^'^Co(s) + c/'Th. (S) ] - c/Pl(5) 



8) 



cTTsCs) ^ 



where Yc = 

Ktm 



TsS-l-i 



<^Tco(5) = | 2 c/'Ts(s)-c/'Th. (s) [l-Kl]} 



These are equations (34) and (35) on page 33. 
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APPENDIX D 



DEVELOPMENT OF THE PRIMARY PIPING TRANSFER FUNCTIONS 



I . Transport Delay Transfer Functions 

The transport' delay of the coolant between the outlet of the reactor 
and the heat exchanger inlet plenum can be expressed as 



where 












' Tho(‘-'^3> 

- Heat exchanger inlet plenum coolant 
temperature (°F) 

' Reactor coolant outlet temperature (°F) 
= Transport time delay (sec) 



(D 1) 



Rearranging terms and expanding equation (D 1) in a Taylor series 






dt 



dt 



(D 2) 



For slow temperature changes, the second order and higher terms 
in equation (D 2) may be ignored. 






+ 13 






dt 



Ho(t) 



(D 3) 



Let 



TH,(t) = T + (t) 

po ip 



(D 4) 



^Ho(^) = ■'h 



(Aho^^^ 



Oo 



where the lowest zero subscript denotes a steady state value and 
the delta prefix a small perturbation about this value. 
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From substitution of equations (D 4) into equation (0 3) and for 



the temperatures in a steady state at t^O 



and 





0 



= c/'Tno(t) = O 



(0 5) 



(0 6 ) 



Impose the perturbations at t=0. Substituting equations (0 4) and 
(D 5) into equation (D 3) with the initial conditions given in equation 
(0 6 ) 

<y’TH,p(t)+ fs ^ (D 7) 

Taking the Laplace transform of equation (0 7) and solving for 

d'lHoO) 1 + fzS 

This is equation (39) on page 36, 

By following a similar development, the transport delay between a 
perturbation in the heat exchanger coolant outlet temperature and the 
resulting perturbation in the reactor inlet plenum coolant temperature 
can be expressed by the transfer function 

(s) _ 1 

(^TcoCs) " l-i-T-lS 

where cTTc-p = Reactor inlet plenum coolant temperature 

perturbation (°F) 

c/'Tco = Heat exchanger coolant outlet temperature 
perturbation (°F) 

= Transport time delay (sec) 
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This is equation (40) on page 



2, Mixing Delay Transfer Function 

Let the reactor and heat exchanger coolant inlet plenums be represent- 
ed by the control volume shown in Figure D-1 . 



Figure D-1. Mixing Volume 

If perfect mixing is assumed, the coolant outlet temperature is equal 
to the average coolant temperature in the plenum. At heat balance on 
the volume requires that the net heat flowing into the plenum be equal 
to the increase in stored energy. Assuming no ambient heat losses, and 
assuming the temperature of the plenum structural material remains 
constant, the heat balance per unit time over a time interval At is 




M = Mass of coolant in plenum (Ibm/sec) 



T = Average plenum coolant temperature 
(°F) 



T. = Plenum inlet coolant temperature 

* / O I- \ 



T^ = Plenum outlet coolant temperature 



m = Coolant mass flow rate (Ibm/sec) 




where C = specific heat of coolant 



The resulting equation is 




(D 8) 
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where T 
Let 



M 

— = mixing time delay, 

m 



T.(t) = T.^ + /T.(t) 
To(t) = To^ t /r„(t) 



(0 9) 



where the lowest zero subscript denotes an initial steady state value 
and the delta prefix a small perturbation about this value. 

Substituting equations (D 8) into equation (D 9) 

-f ^ J'Todt) + </Vo(e) = (DIO) 

Taking the Laplace transform of equation (DIO) and solving for 

J'^^/ cTt, 

O^T.(s) ■ 1+Ts 

This is equation (41) on page 37. 
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INTERACTIVE PROGRAM LISTING 
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LABEL SIMJLATI3N 0= A PRESSURIZED WATER REAC'^OR 
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OUTPUT TIME,XE135,TAy/(500. ,550. ) 

LABEL XE AMD A/3 TEMPERATURE RESPONSE 

LABEL XENON TRANSIENT 

PAGE HEIGHT = 5, WIDTH = 7 
PAGE XVPLOT 
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INTERFACE PROGRAM LISTING 
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